Transient vortex events in the initial value problem for turbulence 
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A vorticity surge event that could be a paradigm for a wide 
class of bursting events in turbulence is studied to examine 
the role it plays in how the energy cascade is established. 
The identification of a new coherent mechanism is suggested 
by the discovery of locally transverse vortex configurations 
that are intrinsically helical. These appear simultaneously 
with strong, transient oscillations in the helicity wavenumber 
co-spectrum. At no time are non-helical, anti-parallel vortic- 
ity elements observed. The new mechanism complements the 
traditional expectation that the development of a peak of the 
maximum vorticity ||w||oo(£) would be connected to nearly 
simultaneous growth of the dissipation, eventually leading 
to the formation of the energy cascade with signatures such 
as spectra approaching -5/3 and strongly Beltramized vortex 
tubes. Comparing how different large-eddy simulations treat 
the spectral transport of helicity demonstrates that the dy- 
namics leading to the helical vortex configurations requires 
both nonlinear transport and dissipation. This finding em- 
phasizes the importance of properly modeling both nonlinear 
transport and dissipation in large-eddy simulations. 
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FIG. 1. Comparison of the growth of 



between our 



256 direct numerical simulation and three results on 64 
meshes. The full 256 3 DNS filtered onto a 64 3 mesh, la- 
beled 64 3 f-DNS, a 64 3 traditional Smagorinsky eddy viscosity 
model, and a calculation using the LANS-a model. This fig- 
ure sets the time scales for our analysis. 

Although three-dimensional turbulence is character- 
ized by intermittent events both in space and time, it 
is often envisaged as a homogeneous, statistically-steady 
tangle of vortex tubes accompanied by a steady trans- 
fer of energy through the spectrum to high wavenumbers 
and a fc -5 / 3 energy spectrum. This statistically-steady 



description has been used to verify turbulence models 
through simulations of forced turbulence and decaying 
homogeneous, isotropic turbulence [QJ|]. In these stud- 
ies, it is sometimes argued that the initial nonequilib- 
rium transients can be ignored as being non-universal. 
However, each intermittent burst of turbulence is itself 
a transient dynamical process involving individual vor- 
tex interactions whose ensemble average is responsible for 
the overall statistical properties. Our objective is to ap- 
ply a coordinated set of diagnostics in both physical and 
wavenumber space for numerically detecting individual 
vortex surge events and their effects upon the subsequent 
dynamics. 

The numerical investigation of the decay of turbulence 
discussed here is designed to improve our understanding 
of individual intermittent events and continues a long 
tradition in numerical modeling of the initial value prob- 
lems in computational fluid dynamics. Numerically, the 
question of transient phenomena in turbulence has been 
considered previously using a variety of initial conditions 
[^-Q. In this letter the evolution from smooth, random, 
initial conditions introduced in H to steady turbulent 
decay is considered using this coordinated set of diag- 
nostics. These diagnostics show that this initial value 
problem for turbulence evolves through several complex 
states in a sequence of transitions. These include: 

• formation of vortex sheets that interact, encounter 
each other transversely and then begin to roll up 
into vortex tubes, 

• development of a peak in the maximum vorticity 
IM|oo(i) that is correlated with helicity signatures 
in both physical and wavenumber space, 

• rearrangement of vortex tubes into transverse pairs 
having oppositely signed helicity A = u • u>, where 
vorticity lo = curlu. Unlike inviscid configurations 
, the velocity u on each tube of the pair is only 
partially induced by its partner. Instead, u on the 
tubes arises primarily as a response to strains and 
dissipation distributed within the pair. 

• formation of the classical decay regime with a k~ 5 ^ 3 

energy spectrum. 
While our analysis also includes traditional diagnos- 
tics such as energy decay, our new understanding will 
arise primarily through: the time evolution of the maxi- 
mum vorticity in Figure [l| physical space visualizations 
in Figures || and ^ helicity probability distributions in 
Figure ^ and the accompanying helicity co-spectrum 
in Figure [| Central to our new understanding is the 
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following observation: although this initial condition is 
not a Beltrami flow, spatial regions develop early dur- 
ing its evolution where the helical alignment measured 
by cos 6* = A/(uu) of either sign is locally near unity in 
magnitude. Evidence is given that the ensuing dynamics 
is influenced by helicity A, whose evolution depends upon 
the interplay between nonlinear transport terms, viscous 
effects and large-eddy stress parameterizations. 

The time scale for this investigation is set by a surge 
in the growth of the maximum vorticity at t — 0.5 in Fig- 
ure [l| This time scale also appears in three 64 3 results: 
the 256 3 DNS filtered onto a 64 3 mesh, in a Smagorin- 
sky calculation using a traditional eddy viscosity model, 
and in a calculation using the new LANS-a model, which 
preserves nonlinear transport properties (see j^] and ref- 
erences therein) . The growth of vorticity comes from the 
vortex stretching terms, which at early times involves 
only large-scale strain and, thus, should not be affected 
by the small scales. Evidence that vortex stretching at 
early times is not strongly affected by the small scales 
is that all three calculations have the same time scale 
t = 0.5 for the vorticity surge. The differences among 
the calculations should therefore be due to how vortic- 
ity is suppressed by either viscous effects or LES model 
effects. The identical growth in \\u\\oo in each case until 
t = 0.3 tells us that until that time dissipation and LES 
parameterization have not yet affected the calculations. 

It has previously been shown that the dominant ini- 
tial structures are vortex sheets that arise out of non- 
Beltrami initial conditions, [||. From the interaction 
of those initial vortex sheets in the weakly dissipative 
regime, Fig. ^ shows that a new configuration of trans- 
versely aligned vortex structures has formed by t = 0.5. 
We will demonstrate that this is an inherently helical 
configuration that has arisen from the non-Beltrami ini- 
tial conditions simultaneously with the vorticity surge 
event. The helical nature of the configuration is shown 
both by the concentrations of helicity density A near the 
vortex structures in Figure [2] and by a skewed, transient 
probability density distribution (PDF) of the cosine of 
the helicity angle at t = 0.5 in Figure ||. Accompany- 
ing this helicity signature in physical space, the helicity 
co-spectrum at low wavenumbers in Figure ^| develops a 
strong signature of alternating sign between wavenum- 
ber bands. This is evidence for the dynamical formation 
of large, asymmetrical distributions of helicity associated 
with the vorticity surge in Fig. |l|. 

While the helicity density is not Galilean invariant, 
Galilean transformations cannot remove the strong fluc- 
tuations in helicity in physical space that we observe, nor 
could they remove the strong fluctuations in the helicity 
co-spectrum in Fig. |[ The implication is that if these 
are configurations that naturally and frequently arise in 
a turbulent flow, then helicity could be playing a central 
role in their dynamical evolution. The PDF of the he- 
licity angle in Fig. || quickly changes into a distribution 
with peaks concentrated at cos 6 — ±1 as reported in 
other turbulent flows [p|Jl^] . The appearance of the in- 



dividual peaks is associated with the formation of nearly 
Beltrami vortex tubes by t — 0.7, shown Fig. ||. 

6j>.55u p] co lines, t=0.5 
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FIG. 2. Isosurfaces of vorticity at t = 0.5 where 
ui > 0.55||lj||oo in red/earth colors with sample vortex lines 
through these regions. Regions of high positive and negative 
helicity are indicated by green and blue respectively. The 
vortex lines meet transversely, which is an inherently helical 
configuration. 
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FIG. 3. Probability density functions of the helicity angle 
cosf? = if ■ a?/(au) in physical space at t = 0.5, during the 
surge in the vorticity maximum, and t = 0.7 after formation 
of the first distinct vortex tubes. The distribution is taken 
over those points with vorticity above the shown threshold, 
55% of the maximum vorticity u) p = ||w||oo, in the subdomain 
containing ||^||oo at both t=0.5 and 0.7. The distribution 
was initially flat. The asymmetry in the distribution emerges 
because the transverse vortex configurations that develop are 
inherently helical. 

The association of the vorticity surge with helicity 
raises some important questions about the absence of 
anti-parallel vorticity elements. In the inviscid limit, 
anti-parallel vorticity elements around the position of 
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the maximum vorticity ||w||ooj were thought to lead to 
the strongest increases in vorticity In a strictly 

anti-parallel configuration there would be equal concen- 
trations of both signs of the hclicity, which would cancel 
at each length scale and thus preclude any spectral oscil- 
lations. 



T = 0.50 




after the vorticity surge at t — 0.5. The helical character 
of these vortex tubes is demonstrated by physical space 
renderings that show that the strongest helicity fluctua- 
tions are clearly tied to vortex tubes in Fig. || at t = 0.7 
and by strong peaks at both 1 and -1 in the helicity dis- 
tribution in Figure ^. This PDF was centered upon the 
domain shown. Note that the two dominant tubes have 
opposite signs of helicity. Rotating one's view of this 
local region shows that the tubes are nearly orthogonal 
and when the entire flow is rendered, this configuration 
appears in a localized corner 1/4 3 of the entire domain. 
This demonstrates that these structures are interacting 
strongly within this localized region. As time goes on, 
vortex tubes whose helicity is opposite develop as or- 
thogonal pairs in many localized regions. This is not the 
first time that nearly orthogonal configurations of vortex 
tubes have been created in direct simulations. They have 
appeared in renderings going back to about mid-1980s 
P,|12| . The new point we are making is that this intrinsi- 
cally helical configuration arises in conjunction with the 
vorticity surge. 



T = 0.70 
.015 1' ' ', ' 




-.015 I I 

2 4 6 8 10 12 
K 

FIG. 4. Comparison of the normalized helicity co-spectra 
A iv (fc) for two times. The normalization is based upon the 
energy and enstrophy at each time so that comparisons in 
magnitudes between the different calculations can be made. 
These strong oscillations first appear at t «s 0.45 when dissi- 
pation is still negligible and the vortex structures are begin- 
ning to change from sheets to tubes via a roll-up instability 
induced by interactions between distinct vortex sheets. The 
higher wavenumber oscillations disappear as soon as dissipa- 
tion becomes significant. 

Instead, at no time during this calculation are the 
the most intense vorticity elements observed to be anti- 
parallel. The calculations that suggest that the anti- 
parallel configuration would be dominant are all based 
upon inviscid calculations using vortex filaments and 
tubes (see jnj and references in Qj). The spectral cal- 
culations reported here, which are viscous and initially 
develop vortex sheets, show no such trend. 

Isolated, helical vortex tubes first form immediately 
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FIG. 5. Isosurfaces of vorticity at t = 0.7 where 
uj > 0.47||cl)||oo in red/earth colors with sample vortex lines 
through these regions. Regions of high positive and negative 
helicity are indicated by green and blue respectively. The vi- 
sualization volume is centered on the new maximum vorticity 
and rotated with respect to Fig. |^ to emphasize Beltrami 
nature of the vortex tubes. From another angle, it can been 
seen that the transverse alignment found in Fig. ^ persists. 

The rapid development of structures characteristic of 
fully-developed turbulence by t — 0.7 from the helical 
state at t = 0.5 occurs simultaneously with the disap- 
pearance of the intermediate peak in the hclicity co- 
spectrum in Fig. ^| seen at t = 0.5. The wavenumber 
(k = 5) of this peak at t = 0.5 is too low for the dis- 
appearance to be a direct effect of viscosity. Instead, 
we believe it is due to a combined effect of transfer 
to larger wavenumbers, followed by dissipation at those 
larger wavenumbers. The evidence for this is the differ- 
ent manner in which the two large-eddy simulations treat 
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the nonlinear transfer of helicity. The LANS-a model 
preserves the transport of helicity. Consequently, when 
energy is transported to small scales and is dissipated, so 
should any helicity that accompanies that energy. The is 
consistent with the observation in Fig. ^ that the inter- 
mediate peak in its helicity co-spectrum disappears by 
t = 0.7. Graphics also show vortex structures consistent 
with the large-scale structures in the DNS. 
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FIG. 6. Comparison of the energy spectra versus the 
three-dimensional wavemimber by shells for several times. 
The main figure shows spectra up to t = 0.9. Between t = 0.3 
and t — 0.5, the high wavenumbers are filled and there is 
little decrease in the overall energy. Between t = 0.5 and 
t = 0.9, there is a significant decrease in the energy in the 
lowest wavenumbers. The expected fc _5//3 law establishes it- 
self by t = 1.3, as shown by the inset. 

For Smagorinsky, the intermediate peak in the helicity 
co-spectrum that had formed by t = 0.5 persists, only 
moving slightly towards higher wavenumbers by t = 0.7. 
There are no assurances in the Smagorinsky model that 
both helicity and energy will be dissipated. In fact, 
Smagorinsky can generate anomalous helicity. Further- 
more, clearly defined vortex tubes have not formed at this 
time in the Smagorinsky calculation. Hence, the proper 
spectral dynamics of helicity and the formation of vortex 
tubes are part of one nonlinear process involving diffu- 
sion, stretching and proper transport. This implies that 
a large-eddy simulation that has correct transport prop- 
erties is required to capture this process. 

Once the interacting, transversely aligned, helical vor- 
tex tubes have formed, dissipation grows rapidly in the 
DNS and as has been reported |5| and shown in Fig. 
the kinetic energy spectrum E{k) gradually approaches 
the classical A; -5 / 3 characteristic of a turbulent energy 
cascade. Once a clear —5/3 spectrum appears near 
t = 1.3, the energy decay becomes self-similar Q. 

Interactions between regions or modes of oppositely 
signed helicity have previously been investigated in the 
context of shell models. In one of the most popular shell 
models, the GOY model, the sign of helicity in each shell 



alternates. It has been shown that the GOY interactions 
are consistent with one channel in a helical decomposi- 
tion of the energy transfer |l^ , |l4| ]. This analysis also 
shows that the strongest wavenumber interactions occur 
between modes of opposite helicity. In DNS calculations 
it has been shown that helicity can arise from non-helical 
initial conditions and that energy and helicity can move 
together to high wavenumbers and be dissipated p^fTsI , 

We conclude that dynamics involving helicity dynam- 
ics is essential during an early stage in establishing the 
turbulent energy cascade from smooth initial conditions. 
The new features that we wish to emphasize are the for- 
mation of a helical state in association with the vorticity 
surge, and how this state rapidly develops into the struc- 
tures and distributions that set the stage for fully de- 
veloped turbulence. We believe that understanding this 
transition could be important in the parameterization of 
transient events in turbulence. Such transitions must in- 
fluence the overall statistical properties controlled by in- 
termittent, intense events. Evidence was presented that 
the mechanism by which classical vortex tubes first form 
requires both transport and dissipation and that this im- 
plies that large-eddy parameterizations should preserve 
some of these helicity transport properties. In ongoing 
work we are investigating the capabilities of several mod- 
ern approaches to large-eddy simulation in representing 
these helicity transport properties. 
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